Anomalous scaling and gapless fermions of d-wave superconductors in magnetic field 
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The problem of d a ,2_ H 2-wave quasiparticles in a weakly disordered Abrikosov vortex lattice is 
studied. Starting with a periodic lattice, the topological structure of the magnetic crystal momenta 
of gapless fermions is found for the particle-hole symmetric case. If in addition the site centered 
inversion symmetry is present, both the location and the number of the gapless fermions can be 
determined using an index theorem. In the case of spatially aperiodic vortex array, Simon and 
Lee scaling is found to be violated due to a quantum anomaly. The electronic density of states 
is found to scale with the root-mean-square vortex displacement as \////(uJ: ms /f 2 ), while thermal 
conductivity is //-independent, but different from the H = case. 



Theoretically, it is well established that the motion of 
the low energy fermionic quasiparticles in the Abrikosov 
vortex state of d— wave superconductors is determined 
by the combination of two effects: the so-called Doppler 
shift [lj and the Z2 Berry phase The former arises 
from the coupling of the quasiparticlc charge current to 
the superfluid velocity and has a classical analog, while 
the latter is purely quantum and corresponds to an extra 
7r Aharonov-Bohm phase accumulated by the quasiparti- 
cle (qp) wavefunction upon encircling a vortex. Semiclas- 
sical analysis of the Doppler shift coupling lead Volovik 
[H to the conclusion that the original nodal spectrum of 
a d-wave superconductor should be filled upon the ap- 
plication of the external magnetic field, H, resulting in 
the finite density of states at zero energy increasing as 
Experimental data on, for example, specific heat 
and thermal conductivity have been largely interpreted 
in this way. 

However, the relatively recent experiments on the ultr- 
aclean YBCO0 have failed to detect the \JtJ dependence 
of the qp contribution to the low temperature thermal 
conductivity n xx . Instead, in the sub Kelvin temperature 
regime k xx /T approaches a value which is essentially H 
independent Q . This is unlike the less pure YBCO sam- 
ples where the said y/H dependence was observed. It 
is worth stressing that the qp mean free path, which has 
been estimated in Ref. 0] , exceeds the magnetic length for 
the H fields above lTesla by several orders of magnitude. 
This suggests that, in the ultraclean samples, the domi- 
nant qp scattering in the Abrikosov state comes from the 
positional vortex disorder, rather than from the impurity 
scattering. This in turn implies, that the quantum effects 
play an important role. 

In this work I suggest a theoretical approach to calcu- 
late the scattering rate due to the weak positional vor- 
tex disorder as well as the resulting thermal conductiv- 
ity. The key idea is to exactly solve the periodic lattice 
problem for a vortex lattice with a site centered inver- 
sion symmetry and a particle-hole symmetry, along the 
lines of Ref. [5(, and then use the symmetry of the re- 
sulting nodal qp wavefunctions to find the effects of var- 



ious perturbations. In this way one can analyze quite 
rigorously the effects of weak positional vortex disorder. 
The new finding is that the coupling of the small dis- 
placement ur of a vortex at R to the new vortex lattice 
nodal quasiparticles [BJ can be thought of as a coupling 
to a small vortex-antivortex dipole. Due to the anoma- 
lous r~2 increase of the exact qp wavefunctions near a 
vortex 0, Q, this coupling is anomalously large and en- 
ters the effective Hamiltonian via the dimensionless ratio 
ur/£ where £ is a short distance cutoff of the order of 
the core radius, rather than via Ur/I as suggested by the 
scaling argument 0, Q . Weakly disordered vortex lattice 
is then viewed as a collection of small vortex-antivortex 
dipoles with random magnitude and random orientation 
whose effect on the new nodal qp can be analyzed within 
the self-consistent Born approximation. The resulting 
vortex lattice disorder-induced density of states is found 
to be finite at zero energy, with an anomalous depen- 
dence on the root-mean-square fluctuation of the vortex 
displacement u rms , scaling as VH f (ul, ns / £ 2 ) , which is 
parametrically larger than in the absence of the anomaly. 
For small u rms , the low temperature thermal conductiv- 
ity is found to be very weakly if-dependent, in qualitative 
agreement with the experimental findings of Ref. 0] . 

The starting point is the lattice Bogoliubov-dc Gcnncs 
(BdG) eigeneauationfioj 



Hoipr = E n ip T . 



(1) 



where the Hamiltonian TIq acts on the two component 
Nambu spinor ip T = [u r ,v r ] T and has the following ex- 
plicit form 
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Both £ T and A r arc defined through their action on a 
lattice function / r as 



4/r = ~t e " lArr+a /r+<5> 
5=±£,±y 
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<5=±i,±y 



In the symmetric gauge, the magnetic flux $ through 
an elementary plaquctte enters the Peierls factor via 
A rr +;r = — 7n/$/</>o, A rr+ £ = 7ra;<f>/0o; the electronic flux 
quantum is (f>o = hc/e. The d-wavc symmetry is encoded 
in r/s = +(— ) if S |j x(y). Orthorhombic distortions can 
be modeled by making horizontal bonds slightly differ- 
ent from the vertical ones, none of which will change the 
key results. The main topological feature of 8 rr > is its 2tt 
winding around the magnetic field induced vortices. The 
initial AnsatzQ for the pair phases is 
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where V x V0(r) = 2ttz £\ 5(r - r ? ) and V • V0(r) = 
where denotes the vortex positions. 

Connecting pairs of vortices by branch cutsj5j, we 
can define the singular gauge transformation^, 0| U = 
e i cr 3<t>r w here the Pauli sigma matrices act on the 
Nambu spinors. The transformed Hamiltonian 7Y(k) = 



-ik-r 



U- 1 H a Ue lk r becomes 

W(k) = a 3 (£ r (k) - n) + <7iA P (k), 



(6) 



where the transformed lattice operators satisfy 



&=±Z,±y 

A r (k)^r = A ^ Z 2,rv+S X TJg^ 6 lfj r+S . (8) 
5=±x,±y 



The physical supcrfluid velocity enters via the factor 



jl _|_ e i(0 r '-0r)| 



(9) 



and represents the lattice analog of the scmiclassical 
(Dopplcr) effect. The Z 2 field Z2,rr' — 1 on each bond ex- 
cept the ones crossing the branch cut where z 2 ,rr' = — 1- 
If the vortices form a regular lattice, the transformed 
Hamiltonian j6]) is invariant under discrete translations 
by the primitive vectors Ri and R 2 defining the mag- 
netic unit cell, reflecting the periodicity of V TT ' and the 
periodic choice of the branch cuts. Consequently, it can 
be diagonalizcd in the Bloch basis. By the Bloch condi- 
tion Tt(k) acts on the periodic functions, and the crystal 
wavevector k varies continuously within the I s * Brillouin 
zone defined by the primitive reciprocal lattice vectors 

K l = 27r z-(RixR 2 ) ' K2 = 2lr i-(R*R 2 ) - 

Note that, since the unitary transformation IA is time 
independent, the Hamiltonians (J5J) and ([6]) have the same 
thermodynamic and tunneling density of states. 

In presence of the particle-hole symmetry (/x = 0) , the 
Hamiltonian which is a In— dimensional Hermitian 
matrix, connects only sites belonging to different sublat- 
tices and can be transformed into the form 



W(k) 



T(k) 
Tt(k) 



(10) 



Following Wen and Zee[ll|,[l2| (see also Ref.QJ]), at all 
points k where the determinant of 7i(k) does not vanish, 
we can define 



MQs) 



H(k) 



= ( ° m \ fill 
|detW(k)|^ U f ( k ) /' ' 



where | det h(k)\ = 1. Therefore, if the phase of dct h(k) 
winds by a nonzero integer multiple of 2n around a loop 
in k-space, then det h(k*) — for some k* inside the loop 
and H(k*) has zero eigenvalues. The associated winding 
number m, can be calculated from 
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where the k-space vector field 



o M (k) = -id klt lndet(/ l (k)) = --Tr [^M^d^M] (13) 

Since a^H* (k)cr2 = — H(— k), we find that the vector 
field is even under k-space inversion, a jU (k) = a p (— k). 
Therefore the circulation of a vanishes at four k points: 
k* = {0,iKi,iK 2 ,i(Ki + K 2 )}. We thus arrive at 
a rather general result, namely if there are zero energy 
eigenstates at any of these four points, they must appear 
in degenerate pairs with opposite m, thereby canceling the 
circulation of a. Such situation is generally unstable to 
cither nucleation of m = +1 and m = — 1 pairs with two 
massless Dirac fermions, or a vacuum of k-space vortices 
implying a spectral gap. 

To proceed with a specific example, consider vortex lat- 
tice with two vortices per magnetic unit cell respecting 
site centered inversion symmetry The space inversion 
symmetry is realized as a projective symmetry and the 
invariant gauge group fl2| is Z 2 , i.e. ordinary space inver- 
sions X must be followed by Z 2 gauge transformations 7 r 
By constructing an operator V = (— l) x+v ~f r T which 
anticommutes with 7i(k) at the four k* points, one can 
show that the spectrum is gapless at these points with 2 
degenerate Dirac cones at each k* giving the total of 8 
massless Dirac fermions in the I s * Brillouin zone [||. 

What happens if the site centered inversion symmetry 
is broken? According to the above general classification, 
cither the spectrum is gapped, or the k-space vortex- 
antivortex pairs are nucleated resulting in two massless 
Dirac fermions. As shown below, the latter case is re- 
alized for weak breaking of the site centered inversion 
symmetry. 

Without loss of generality, let us concentrate on H(0). 
Since H(0)V = -VH(0) and since TrV = 4, the four zero 
energy eigenstates of 7i(0) are also eigenstates of V with 
the same eigenvalue [HJ of +1. It is also straightforward 
to see that two of them can be chosen to be even under 
the projective inversion 7 r X and the remaining two to be 
odd. The nullspace basis is then 



>t : {|1>, |2> = ier a |l>*} - - - etien 
B:{|3),|4) =ia 2 \3)*}...odd, 



(14) 
(15) 



where the first set is even under -y r T and the second one is 
odd. Since these states are +1 eigenstates of V, they are 
also eigenstates of (— l) x+v , with the eigenvalue +1 for 
the even ones and —1 for the odd ones. This means that 
states which belong to the representation A are non-zero 
only on one sublattice, while the states which belong to 
the representation B are non-zero on the complementary 
one. 

The dispersion of the energy near the degenerate 
quadruplet can be found within degenerate perturbation 
theory. To the first order, the low energy effective Hamil- 
tonian is given by the nullspace matrix elements 



7^(kiUi,u 2 ) = (m|W(k;ui,u 2 )| n) , 



(16) 



where we explicitly included the dependence of the 
Hamiltonian on Uj's which denote the possible small dis- 
placements of the positions of the two vortices within 
the magnetic unit cell. Due to its particle hole symme- 
try, 7i(k; Ui, 112) cannot mix states with the same parity 
and therefore Hf££ is block off-diagonal with the struc- 
ture (fTT)|) . By using cr 2 H*(k)cr 2 = —H(—k), one can 



show that Till 0*5 u i, u 2) = — TC 2 \ (~k;ui,u 2 ) and 
(k;ui,u 2 ) = H23^(— k; Ui, u 2 ). In addition, the 
projective inversion 7 r X guarantees that 

Htii(k; u l5 u 2 ) = -H e Jf(-k; -u 2 , - Ul ). (17) 

Clearly, for Uj = we have two degenerate gapless 
Dirac fcrmions at k = and vanishing circulation of a 
(|12| as required by the above argument. 

For finite Uj the zero energy states exist provided that 



detW e// (k;ui,u 2 ) = 



(18) 



Due to (fT7| . we can expand near k = 0, Uj ; = to find 
Til f 3 f (k;u 1 ,u 2 ) = k^A^ +uj fJ ,aj l and H e /J (k; m, u 2 ) = 
k^B^+Uj^. Note that = a 2 = a p and 6* = = 
There are two solutions to the Eq. (fT8"]) , at 



+ \ui a b l a \ 



\A a A l a/3 iiii3\ 2 + \B a K l a(j Uit3 
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where A^„ ~ Ze^x^sm \A\a^* v + B\W*J\ is independent 
of Uj. We see that small vortex displacements induce a 
k-space vortex-antivortex pair, resulting in a pair of gap- 
less fermions near each k*. By going to the second or- 
der degenerate perturbation theory, I have checked that 
the character of the solution does not change. This is 
a result of the general topological structure ofH(k): the 
spectrum is gapless and the number of the Dirac fermions 
is invariant under a finite range of periodic perturbations 
which can break site centered inversion, but do not break 
particle-hole symmetry. In other words, any p-h sym- 
metric perturbation, such as supercurrent variations and 
bond or pair density waves, must be sufficiently strong 



to cause the annihilation of the k-space vortex-antivortex 
pair. 

We now turn to the analysis of the magnetic field de- 
pendence of the effective Hamiltonian (|16|) . According 
to the scaling argument of Simon and Lee 8], since the 
vortex displacement Uj has a dimension of length, in the 
limit of large £ we should haveQ 



1 



^mra (k; Uj ; £) — pTmn' \^k 



£ 



where is the scaling function. As a result the ef- 

fective velocity of the gapless Dirac fcrmions should be 
£ independent, and in addition, the vortex displacement 
fields Uj should enter the effective Hamiltonian as Uj / £ 2 , 
making a?„ ~ V„ ~ £~ 2 . 

Remarkably, explicit numerical evaluation of the ma- 
trix elements ([To]) shows that while the effective velocity 
is indeed ^-independent Q , the vortex displacements cou- 
pling to the gapless fermions is much stronger and instead 



£ . The correct scaling form is rather 



— 1 -p(anom) ( fu. }h_ 

£ mn \ ' Co 



(19) 



where £0 is a short distance cutoff of the order of the 
vortex core radius. This is a consequence of a quantum 



anomaly [lj, [15[ which arises from two effects: first, due 
to the 1/r increase of the superfluid velocity near a vor- 
tex, the extended quasiparticle wavefunctions experience 
an anomalous increase0,[ll,[l3| ~ 1/Vr£ cut off only by 
a short distance scale £o- Second, the potential due to 
the displaced vortex, as experienced by the nodal quasi- 
particle of the unperturbed vortex lattice, is equivalent 
to the potential due to vortex-antivortex " dipole" . At 
distances much larger than |uj|, this "dipole" potential 
falls off as ~ 1/r 2 . The resulting matrix element scales 
as I u j I / (££,0), rather than |u,,-|/£ 2 . The scale invariance of 
the linearized Dirac equation^ is thus broken and just 
as in the theory of critical phenomena, the short distance 
cutoff effects long distance properties via an anomalous 
dimension. In the case at hand, the anomalous dimension 
of Uj is effectively 1 . 

In the weakly disordered vortex lattice, the disorder- 
free vortex lattice nodal quasiparticles experience a ran- 
dom potential produced by the vortex dipoles with ran- 
dom orientation and strength. The vortex displacements 
u Rj > corresponding to the vortex position Rj , can be 
assumed to be uncorrelated 



The effective Hamiltonian near k = is TL^L 



(20) 
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where v„ = 



-B 
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and the 



a M b^ 
L 2 is the area of the 



-l 



f \ R* A* 

sum is over the vortex positions 
system, and just as before, a M ~ 6^ ~ £ 

Near k = 0, we can treat the vortex- antivortex dipole 
scattering within the self-consistent Born approximation 
and the quasiparticle (retarded) self-energy 



E(o;) - 2u 2 rms [\a\ 2 + \b\ 2 ] £■ 



^0(°>( W -E( W ))(22) 



where the integral is over the magnetic Brillouin zone, 
and H 2 + \b\ 2 = E„(M 2 + K\ 2 ); G { q } is the quasipar- 
ticle Greens function for the disorder free vortex lattice. 

We find the vortex disorder scattering rate, r" 1 , from 
the imaginary part of the retarded self-energy to be 
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where uy j_ are the velocities of the dispersion of the 
(disorder-free) vortex lattice Dirac quasiparticles along 
the principal axes. Note, that due to the quantum 
anomaly, the scattering rate has the scaling form 



'I 2™a 
so 



This should be contrasted with the result obtained from 
the naive scaling argument T/jn = VHT (uf ms H/$>o) , 
which is parametrically smaller than (|23|) . 

The contribution to the zero energy density of states 
from the vicinity of k = is therefore 



A^o(0) = 



In (TTT vy /VfIIi x ) 



TT V±V\\T l 

Upon summing the contribution from all four nodes, 
we find that that the vortex disorder induced density of 
states has a scaling form iV(0) = V~Hp(u 2 ms /t; 2 ). 

Since the vortex disorder induced self-energy was ar- 
gued to be given by the self-consistent Born approxima- 
tion, in order to find the electronic contribution to the 
thermal conductivity k xx we need to calculate the vertex 
correction by summing the ladder diagrams. In the limit 
of small u rms we find 
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(24) 



where (3 V = 1 + 0{u* ms ) > 1. 

The above result is independent of magnetic field and 
depends only on the vortex lattice renormalized Dirac 
cone anisotropy. The extra factor of 2 relative to Ref . [IH 
comes from the node doubling Q ■ 

In conclusion, I showed that weak vortex lattice disor- 
der exhibits an anomalously strong coupling to the low 
energy fcrmions, which violates Simon-Lee scaling in a 



novel way. The resulting vortex disorder induced density 
of states was calculated for weak disorder and found to 
scale as \/Hf(u 2 ms /£ 2 ). The electronic contribution to 
k xx /T was found to be magnetic field independent, but 
different from its H = value. The latter can be under- 
stood to result from the cancelation of the 71-dependence 
of the vortex disorder induced density of states and the 
mean free path, which is set by the intervortex sepa- 
ration in this ultraclean limit. This result is in quali- 
tative agreement with the experimental finding on the 
ultraclean YBCOQ, although quantitative agreement is 
unlikely due to the significant vortex disorder in the ac- 
tual sample. Nevertheless, the theory presented here, 
while strictly limited to the case of weak positional vor- 
tex disorder, does suggest that as u rms increases and 
becomes of order of the magnetic length £, the vortex 
disorder induced density of states continues to scale as 
s/H. Since the mean free path is still set by £, this sug- 
gests TJ-independcnt, albeit vortex disorder dependent, 
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